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The high statistics data from the pp — > ppr\ reaction measurement, delivered by the COSY- 
1 1 collaboration, are now being evaluated using the correlation femtoscopy technique. This 
method is based on the relative momentum correlations of two emitted protons and may per- 
mit determination of the size of the reaction volume. 

For the very first time, we apply an intensity interferometry technique to study the mechanism 
of the meson production via the nucleon-nucleon interaction close to the kinematical thresh- 
old. We invented a method to determine correlation function for the ppr\ system free from the 
physical multi-pion production background. We show the comparison of experimental results 
with theoretical predictions and appraise the accuracy achieved for the determination of the 
size of the emission source. 

PACS: 13.60.Le, 14.40.-n, 25.75.Gz 

1 Introduction 

One of the possible mechanisms of the rj meson production via the pp — > ppr\ reaction is a di- 
rect emission of the r\ meson from the interaction region and the other possibility, believed to 
be dominant, is the creation of this meson via the resonant state S'n(1535) |Q]|2E). Very close 
to the kinematical threshold the direct influence of the resonant state on the form of differential 
distribution of the cross section (e.g. Dalitz plot) is very difficult to observe due to the fact that 
the resonance is more than an order of magnitude broader (r« 150 MeV |4|) than the available 
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range of the invariant mass distribution, which in the case considered in this article amounts to 
15.5 MeV 0. Therefore it is worth trying to find an observable different from the distributions 
of the cross sections which can deepen our understanding of the reaction dynamics. We assess 
that such an additional information about the production mechanism may be derived from the 
size of the reaction region from which the particles emerge. 

In comparison to the direct production, the presence of the resonance in the production pro- 
cess (pp — > pN* — > ppij) will increase the effective size of the protons separation at the moment 
when they appear as free (on-shell) particles. The information about the effective spatial pro- 
ton's separation at this moment may be derived from the proton-proton correlation function. 
However, the shape of this function depends also on the interaction between the emitted parti- 
cles, and in order to draw unbiased conclusions about the spatial size of the emission source, 
the effect from the final state interaction (FSI) has to be determined precisely. Due to the pres- 
ence of the Coulomb and strong interaction, and the non-trival problems with the description 
of the three body (nucleon-nucleon-meson) system, there exists at present no rigorous descrip- 
tion of the interaction within the ppr\ system Yet, in practice we can account for the ppr\ 
FSI, parametrizing it from the precisely measured m m and m pp invariant mass distributions 1 7 1 . 
Therefore, by the analogy to the successful analysis based on the Dalitz plot and correlation func- 
tion, performed for the determination of the size of the neutron halo in the 14 Be ( 12 Be + n + n) 
and the mechanism of its dissociation on C and Pb targets 1 8 1, we would like to conduct similar 
investigations for the pp-q system in pp — > ppr\ reaction. In case of the pp — > pp-q reaction it can 
be possible that two protons coming from the reaction volume could be emitted with a uncorre- 
lated phases since the reaction proceeds in two steps (pp — » pN* — > pprf) via the excitation of 
the intermediate resonance state N* (1535) which decays delayed in time into the proton and the 
r] meson. 

Due to the short life-time of the resonant state, the average increase of the effective size of the 
source in respect to the direct production is expected to be small (in the order of the fraction of 
femtometer), but as we will argue later, the accuracy we can achieve is better than 0.1 fm. More- 
over, there is another interesting aspect connected to the study of the reaction volume. Namely, 
it cannot be a priori excluded that the geometrical dimensions of the source can be much larger 
due to the reflection of the topological (Borromean) bounding in the pp-q system. As Borromean 
we call a bound three-body system in which none of the two-body subsystems is bound. There 
exists, an encouraging and interesting conjecture of S. Wycech who pointed out that a large en- 
hancement in the excitation function of the pp — > pp-q reaction observed close to the kinematical 
threshold may be described assuming, that the proton-proton pair is produced from a large (Bor- 
romean like) object of a 4 fm radius 1 9 1 . 

The Borromean bounding in the three body system is a very intriguing issue attracting the 
researchers from different fields of science. In nuclear physics the 11 Li ( 9 Li + n + n) and 6 He 
(a + n + n) nuclei have been found to have such a property 1 10 1. Another, interesting and worth 
noticing example, is that very recently nanoscale Borromean rings were constructed in a wholly 
synthetic molecular form II 1111 21 . 

At present it is, however, still not established whether the low energy pp-q system can form a 
Borromean or resonant state. 

In the subsequent sections we will present a correlation function for the ppr\ system produced 
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in the collisions of nucleons, where the r\ meson is identified via the missing mass technique. In 
such measurements it is impossible to identify the production of ppr] system on the event-by- 
event basis due to the physical background originating from the pp — > pppions reactions. We 
have, however, succeeded to develop a method for constructing a correlation function free from 
the multi-pion production background. The technique will be described in this article. 



2 Correlation femtoscopy 

The momentum correlations of particles at small relative velocities are widely used to study 
space-time characteristics of production processes, so serving as a correlation femtoscopy 1 13 1. 
Correlation femtoscopy originates from an intensity interferometry known as HBT effect 1 14 1. 
Implemented into nuclear physics I15I16I13I . this technique permits to determine the size of the 
source from which the particles are emitted. It is based on the correlation function, which relates 
the space-time separation of the particles to their momenta p\ andj<2 at the emission time. This 
function can be expressed in terms of pair- and single-particle cross section 1 17 1: 

R{P1 ' P2) = ° ■ (d^/d^^a/d^) ~ l! (1) 

where C denotes the overall normalization constant. Generally one can relate a correlation func- 
tion with a Fourier transform of the spatial distribution of the emission source 1 18 1. The shape 
of a two-proton correlation function depends on the spatial size of the source, quantum statistics 
and the interaction between the outgoing particles. 

In our case the influence of the ppr\ FSI can be derived from the experiment. More details will 
be given in the next section. Thus, taking FSI into account when calculating a correlation function 
will permit to study the dependence of the shape of this function on the spatial dimensions of the 
reaction region. 

As far as the experiment is concerned, in case of impossibility for the single particle cross 
section measurement, an alternative definition for the correlation function can be applied. Equiv- 
alently to the equation^ one can calculate such a function by generating an "uncorrelated" yield 
via event mixing technique, as shown by Kopylov and Podgoretsky 1 16 1. The correlation function 
is then given by 1 17 1: 

where Y\2{q) denotes the coincidence yield and Y* 2 (q) stands for the yield derived from the 
uncorrelated reference sample. C* 2 is an appropriate normalization constant. Here, R(q) denotes 
a projection of the correlation function onto the relative momentum of emitted particles q = 

\fi -f%\- 



3 Interaction within the ppr] system 



The influence of the interaction between ejectiles of the pp — » ppr/ reaction on the correlation 
function can be inferred from the invariant mass distributions of the two-particles subsystems. 
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Fig. 1. Missing mass spectrum for pp — > ppX rection determined in the experiment at beam momentum 
Pb = 2.0259 GeV/c |5| (points). The solid line histograms present the simulation of 1.5 • 10 8 events of 
pp — > ppr\ reaction, and 10 10 events for the pp — + pp2n, pp — » pp37r and pp — > pp47r reactions. The 
simulated histograms were fitted to the data varying only the magnitude. 



The pp — > ppr/ reaction was measured with high statistics by the COSY- 11 collaboration at 
beam momentum ps = 2.0259 GeV/c |5|. The determined missing mass spectrum for the 
pp — > ppX process is presented in figure ^ One can easly recognize the sharp peak from the 
r\ meson production over a flat multi-pion production background. The estimated number of r\ 
production events is around 25 thousands. High statistics in an such experiment allows to derive 
the distributions of differential cross sections free of the multi-meson production background 
J5][7). An example of such a distribution is presented in the left panel of figure [2] This figure 
presents the projection of the phase-space distribution onto the two-proton invariant mass axis. 
One can easily recognize the growth of the population density at the range where protons have 
small relative momenta. The shown distribution is free of the multi-pion background since the 
number of pp — > ppr\ events has been extracted for each point separately, first selecting data 
according to the square of the proton-proton invariant mass (s pp ) value and then for each s pp 
interval constructing a missing mass distributions. Generally, the full experimentally accessible 
information of the interaction within the ppr\ system is contained in the Dalitz plot distribution 
(shown in the right panel of figure|2j which we can use when taking into account an influence of 
the ppi] FSI on the correlation function. 

4 Proton-proton correlation function for the pp — > pprj reaction 

In this section we would like to present a method of evaluation of proton-proton correlation 
function for the investigated pp — * ppr\ reaction. Subtraction of the physical background and 
corrections for the limited acceptance of the detection system constitute the two main challanges 



Correlation femtoscopy for studying i] meson production mechanism 



5 




Fig. 2. (Left panel): Distribution of the square of the proton-proton (s pp ) invariant mass for the pp — » ppr\ 
rection at Q — 15.5 MeV 15 171 . The solid line corresponds to the calculations under the assumption of only 
proton-proton interaction 1 19 1. (Right panel): Dalitz plot distribution. Experimental results determined for 
the pp — > ppr\ rection at the excess energy of Q = 15.5 MeV. Data were corrected for the detection 
acceptance and efficiency. 

we have to face when deriving the correlation function from the experimental data. Here we will 
concentrate on the former one. 

According to the equation|2]the two-proton correlation function R(q) can be defined 11171 as 
a ratio of the reaction yield Y ppv (q) to the uncorrected yield Y*(q). In the discussed pp — > ppij 
experiment, only four-momenta of two protons were measured and the unobserved meson was 
identified via the missing mass technique 115 17 1 . In such measurement the entire accessible infor- 
mation about an event is contained in the momentum vectors of registered protons. Therefore, it 
is in principle impossible to decide whether a given event corresponds to the rj meson production 
or whether it is due to the multi-pion creaction. However, statistically, on the basis of the missing 
mass spectra, one can derive a number of events originated from the production of ppr\ system, 
for a chosen region of the phase-space. Therefore, Y ppv (q) can be easily extracted for each stud- 
ied interval of q by dividing the sample of measured events according to the value of q, next 
calculating the missing mass spectra of the pp — > ppX reaction for each sub-sample separately, 
and counting the number of ppr\ events from these spectra. An example of such histogram for 
one value of q is presented in the right panel of figure[3] To demonstrate the full set of the data 
in the left panel of figure [5] we show also a spectrum of the population density over the plane 
spanned by q and the missing mass. 

The statistics obtained in the considered measurement allowed to divide the kinematically 
available range of q into bins whose width (Aq = 5 MeV/c) corresponds approximately to the 
accuracy of the determination of the relative proton momentum (cr(q) ~ 6 MeV/c). 

An extraction of Y*(q), unbiased by the multi-pion production is however not trivial. Ap- 
plying a mixing technique one can construct an uncorrected reference sample taking momentum 
vectors of protons corresponding to different real events. A real event is determined by the mo- 
mentum vectors of two protons registered in coincidence, and an uncorrelated event will thus 
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Fig. 3. (Left panel): The example of the high statistics distribution of the measured pp — > ppX reaction 
close to the r\ meson production threshold |5|. (Right panel): The example of missing mass spectrum 
measured 1 5J for the pp —> ppr\ reaction at q 6 (0.06; 0.07) GeV/c. 



comprise momentum vectors of protons ejected from different reactions. Unfortunately, in such 
a sample of uncorrected momentum vectors, due to the loss of the kinematical bounds, the pro- 
duction of the T] meson will be not reflected on the missing mass spectrum and hence it can not 
be used to extract a number of mixed-events corresponding to the production of the r\ meson. 
Therefore, in order to determine a background-free correlation function for the pp — > ppr\ reac- 
tion we performed an analysis in the following manner: First, for each event, we have determined 
a probability uj that this event corresponds to the pp — > pprj reaction. The probability w,, that 
i th pp — ► ppX event with a missing mass TOj, and relative momentum of qi corresponds to 
pp ~ * ppr\ reaction was estimated according to the below formula: 

N 

ojj = — fm,-,a,), (3) 

N V + N B ,q h y ' 

where N„ stands for the number of the pp — > ppr\ reactions and Ng is the number of events 
corresponding to the multi-pion production. The values of N v and Nb were extracted from 
the missing mass distributions produced separately for each of the studied intervals of relative 
momentum q. An example of a missing mass spectrum with the pictorial definitions of Nb and 
N v is presented in the right panel of figure [5] Next, in order to obtain a value of Y pprj (q) for a 
certain q we have added probabilities of all events for which a relative momentum of two protons 
belongs to the bin centered at q. Thus, Y ppri (q) = ^ uji where i enumerates all events with 
relative protons momentum corresponding to the considered q bin. The result is of course per 
definition the same as obtained earlier by counting the number of events (Y pprj (q) = J2j N v (j)) 
where j enumerates bins in the histogram. 

However, now having introduced the weights we can calculate also the value of Y* (q) with- 
out a bias of the multi-pion background. We can achieve this by sorting an uncorrelated sample 
according to the q values similarly as in the case of the correlated events and next for each sub- 
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Fig. 4. Comparison of the experimental correlation function for the pp — > ppX reaction, without ac- 
ceptance correction represented by full dots and theoretical calculations indicated as two solid lines for 
emission sources parametrized by Gaussian distribution with r = 2.0 and 4.0 fm. 



sample we construct background free Y* (q) as a sum of the probabilities that both protons in an 
uncorrected event originate from the reaction where the rj meson was created. Specifically, if in 
a given uncorrected event denoted by k, one momentum is taken from a real event say kl and the 
second momentum from an real event k2, then the probability that both correspond to the reac- 
tions where the r\ was created equals to ojki ■ ^k2, and hence the uncorrected yield Y * (q) may be 
constructed as a J2k Ukl ' Wfe 2, where k enumerates events in the uncorrected sub-sample with 
a relative protons momentum q. 

In the limit of absence of background all values of u> would be equal to 1 and so deriving 
both Y ppv (q) and Y*(q) we would just count the events. In the presence of the background if 
a missing mass of real event is from the outside of the -q signal than it is certain that this event 
corresponds to the multi-pion production process and it indeed will not be taken into account 
in the calculation of Y ppv (q) and Y*(q) since its weight will be equal to zero. If the missing 
mass of the given event is from the region where N v is larger than zero than we don't know 
whether the event is from the r\ or multi-pion creation. However, for the calculation of the yields 
we need only the overall number of i] events from a given range of the q value. This could be 
determined either by taking only r\ events into account and counting them with a weight equal to 
one (which is in principle impossible since this events are kinematically indistinguishable from 
the multi-pion production) or as we did by counting all events fulfilling a required kinematical 
conditions (in our case we look only at q) and assigning them probabilities that they belong to 
the r] production. 
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The correlation function derived from the data is presented in figure 0] It is compared to the 
calculations, performed assuming a simultaneous emission of the two protons and the rj meson 
and approximating tentatively the effective spatial shape of the emission zone by the Gaussian 
distribution. In such a case the standard deviation of such distribution - hereafter referred to as 
ro - constitutes the measure of the dimension of the source. For the simulations we presently 
adapted a computing procedure written by R. Lednicky 12012 II . which already has been success- 
fully applied in other studies (example in ref. 1 8 1). At present we take into account proton-proton 
FSI only, but in further analysis, we would like to make simulations with compliance of the pprj 
FSI as can be extracted from the Dalitz plot distributions. 

The determined experimental correlation function shows a maximum at a value of q as pre- 
dicted by simulations. The height of the peak at q w 40 MeV/c depends significantly on the 
value of ro and therefore the magnitude of this maximum may serve as a measure of the volume 
of the reaction zone. A rough comparison between theoretical correlation functions simulated 
with ro = 2 fm and 4 fm, respectively, and the experimental one shown in the figure |4] indi- 
cates that the size of the reaction volume can be approximated by the Gaussian distribution with 
ro w 3 fm. Of course, generally the size of the reaction zone can be extracted by the comparison 
of the experimental and simulated correlation functions treating tq as a fitting parameter. 

The results presented in the figure^have not been corrected for the acceptance of the COSY- 
1 1 detection setup yet, and the theoretical calculations have been performed without taking into 
account the experimental spread of the momenta. These corrections could significantly influence 
the final results. Therefore in a present article we do not perform a quantitative analysis. We 
only put out the emphasis on exploring the method which can be further adapted to examine the 
dynamics of rj and rf meson production mechanism in nucleon-nucleon collisions. Although 
after the mentioned corrections the result may be different, it is worth noticing that the extracted 
value of ro is compatible with the effective range (2.8 fm 1221 'I of the proton-proton interaction. 

From the comparison of the theoretical lines and the data in figure |4] one can appraise a 
statistical accuracy of the ro determination which can be achieved when performing fit with ro as 
a free parameter. Examining in figure0]the range of 0.02 GeV/c < q < 0.10 GeV/c one observes 
there 15 points with a statistical errors corresponding to the accuracy varying between 0.2 fm 
and 0.5 fm. Thus as far as the statistical precision is concerned we will be able to determine the 
effective size of the source with the precision much better than the errors quoted above for the 
single points. 

Acknowledgement: We acknowledge the support of the European Community-Research Infras- 
tructure Activity under the FP6 "Structuring the European Research Area" programme (Hadron- 
Physics, contract number RII3-CT-2004-506078), of the FFE grants (41266606 and 41266654) 
from the Research Centre Jiilich, of the DAAD Exchange Programme (PPP-Polen), of the Polish 
State Committe for Scientific Research (grant No. PB1060/P03/2004/26), 
and of the RII3/CT/2004/506078 - Hadron Physics-Activity -N4:EtaMesonNet. 



Correlation femtoscopy for studying i] meson production mechanism 



9 



References 

[1] P. Moskal, M. Wolke, A. Khoukaz, W. Oelert: Prog. Part. Nucl. Phys. 49 (2002) 1. 

[2] G. Faldt, T. Johansson, C. Wilkin: Phys. Scripta T 99 (2002) 146. 

[3] C. Hanhart: Phys. Rept. 397 (2004) 155. 

[4] S. Eidelman et al.: Phys. Lett. B 592 (2004) 1. 

[5] P. Moskal et al.: Phys. Rev. C 69 (2004) 025203. 

[6] F. Kleefeld: Schriften des FZ-Jiilich, Matter and Material 1 1 (2002) 51, 

e-Print Archive: nucl-th/0 108064 , and this proceedings, e-Print Archive: nucl-th/05 10017 

[7] P. Moskal: e-Print Archive: hep-ph/0408162 

[8] F. M. Marques et al.: Phys. Rev. C 64 (2001) 061301(R). 

[9] S. Wycech: Acta Phys. Pol. B 27 (1996) 2981. 

[10] M. V. Zhukov et al.: Phys. Rept. 231 (1993) 151. 

[11] K. S. Chichak et al.: Science 304 (2004) 1308. 

[12] S. J. Cantrill et al.: Acc. Chem. Res. 38 (2005) 1. 

[13] R. Lednicky: Nukleonika 49 (Sup. 2) (2004) S3. 

[14] R. Hanbury-Brown, R. G. Twiss: Phil. Mag. 45 (1954) 663. 

[15] S. E. Koonin: Phys. Lett. B 70 (1977) 43. 

[16] G. I. Kopylov, M. I. Podgoretsky: Sov. J. Nucl. Phys. 15 (1972) 219. 

[17] D. H. Boal et al.: Rev. Mod. Phys. 62 (1990) 553. 

[18] G. Baym: Acta Phys. Pol. B 29 (1998) 1839. 

[19] V. Baru et al.: Phys. Rev. C 67 (2003) 024002. 

[20] R. Lednicky: private communication (2005). 

[21] R. Lednicky and L. Lyuboshits: Sov. J. Nucl. Phys. 35 (1982) 770. 

[22] J. P. Naisse: Nucl. Phys. A 278 (1977) 506. 



